Next-leading BFKL effects in forward-jet production at HERA 
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We show that next-leading logarithmic (NLL) Balitsky-Fadin-Kuraev-Lipatov (BFKL) effects 
can be tested by the forward-jet cross sections recently measured at HERA. For da/dx, the NLL 
corrections are small which confirms the stability of the BFKL description. The triple differential 
cross section da / dxdk^dQ^ is sensitive to NLL effects and opens the way for an experimental test 
of the full BFKL theoretical framework at NLL accuracy. 



1. Forward-jet production in deep inelastic scatter- 
ing is a process in which a jet is detected at forward 
rapiditites in the direction of the proton. The virtuality 
of the intermediate photon and the squared trans- 
verse momentum of the jet are two hard scales. When 
the total energy of the photon-proton collision W is as- 
sumed to be large, corresponding to a small value of the 
Bjorken variable x, forward-jet production is relevant jlj 
for testing the Balitsky-Fadin-Kuraev-Lipatov (BFKL) 
approach [2]. 

Indeed in the small— a: regime, besides the large log- 
arithms coming from the strong ordering between the 
proton scale and the forward-jet scale (which are re- 
summed using the Dokshitzer-Gribov-Lipatov-Altarelli- 
Parisi (DGLAP) evolution equation Q), other large log- 
arithms arise in the hard cross section itself, due to the or- 
dering between the energy W and the hard scales. These 
can be resummed using the BFKL equation, at leading 
(LL) and next-leading (NLL) logarithmic accuracy 
By contrast, in fixed-order perturbative QCD calcula- 
tions the hard cross section is computed at fixed order 
with respect to as, and the next-to-leading order (NLO) 
predictions fail to describe the data. 

NLL corrections to the LL-BFKL kernel were found to 
feature spurious singularities. It has been realised that 
renormalisation-group improved NLL regularisations can 
solve this singularity problem and lead to consistent 
NLL-BFKL kernels. This, along with the success 0, ll| 
of the LL-BFKL approach in describing the da/dx data, 
motivates the present phenomenological analysis of NLL- 
BFKL effects in forward-jet production. 

When fitting da/dx, we obtain that the NLL correc- 
tions are small, which results in a good description of 
the HI data by NLL-BFKL predictions. We also show 
that the recently measured triple differential cross sec- 
tion da/dxdk^dQ^ ^ allows for a detailed study of the 
NLL-BFKL approach and of the QCD dynamics of for- 
ward jets. In particular, it has the potential to address 
the question of the remaining ambiguity corresponding 
to the dependence on the specific regularisation scheme 
of the NLL kernel. 



The present study is a phenomenological analysis of 
the new forward-jet data using NLL-BFKL predictions 
depending on the regularisation schemes and the renor- 
malisation scale. In Ref. [13], such a phenomenological 
investigation has been devoted to the proton structure 
function data, taking into account NLL-BFKL effects 
through an "effective kernel" (introduced in ^]) using 
three different schemes. A saddle-point approximation 
for hard enough scales is used in order to obtain a phe- 
nomenological description of NLL-BFKL effects. In the 
present study devoted to forward-jet production, we im- 
plement them in a similar way. This allows to study the 
NLL-BFKL framework even though the determination of 
the next-leading impact factors is still in progress 

2. We shall use the usual kinematic variables of deep 
inelastic scattering: x = Q'^ / {Q'^ -\-W^) and y — Q'^/{xs) 
where ^/s is the center-of-mass energy of the lepton- 
proton collision. In addition, xj is the jet longitudinal 
momentum fraction with respect to the proton. The fully 
differential cross section for forward-jet production reads 



dxdQ'^dxjdki, 



hiy) 



dal*^-"" 
dxjdk^ 
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with friy) = l-y + y'^/2, /^(y) = 1 - y, and where 
da^*f/~*'^'^ /dxjdk^ is the cross section for forward-jet pro- 
duction in the collision of transversely (T) or longitudi- 
nally (L) polarized virtual photons with the proton. 

In the following, we consider the high-energy regime 
X <C 1 in which the rapidity interval Y = \og{xj/x) is 
assumed to be large. Following the phenomenological 
NLL-BFKL analysis of 0, one obtains for the forward- 
jet cross section: 



da. 



T,L 

dxjdk^ 



as{k'^)as{Q'^) 



feff{xj,kl) 



c?7 



0T:,^(7) e"('=^'5)x.//(7,a(fcTQ))i' (2) 



with the running coupling constant given by 

a(fc2) = as{k')Nj7r = [6 log (A:VA|ci?)] " 



(3) 



2 



- ^ n 

0.8 - 
0.6 



CCS 

S3 

S4 



LL BFKL 




0.05 0.1 



FIG. 1: 7c, QXe//(7c,a) and aXe'//(7c, a), as functions of 
a for the three BFKL resummation schemes CCS, S3 and 
S4. The fixed LL values are respectively 1/2, 0.43, and 5.47 
corresponding to a = 0.16. 



where 6= ll/12-7V//6iVc. 

In formula ([2]), feff{xj,k^) = g+{q + q)CF/Nc is the 
effective parton density which resums the leading log- 
arithms in log(/c^/AQ(^^). g (resp. q, q) is the gluon 
(resp. quark, antiquark) distribution function in the inci- 
dent proton. Since the forward jet measurement involves 
large values of kx and moderate values of xj, formula 
([2]) features the collinear factorization of feff, with k'^ 
chosen as the factorization scale. 

The NLL-BFKL effects are phenomenologically taken 
into account by the effective kernel Xe//(7,Q:). Indeed, 
the NLL-BFKL kernels Xnll{7,^) provided by the reg- 
ularisation procedure obey a consistency condition [I4I 
which allows to reformulate the problem in terms of 
Xeffil, a)- We shall consider the CCS scheme [f|] and the 
S3 and S4 schemes [Hj in which xnll also depends ex- 
plicitly on a. In each case, the effective kernel Xeffil, 01) 
is obtained from the NLL kernel Xnll{,1,^) by solving 
the implicit equation 



Xeff = XNLLil^a Xeff) 



(4) 



as a realization of the consistency condition. 

The details of the approximation ^ are given in [lo| . 
with the only difference that the kernel considered for F2 
was naturally the asymmetric one xjvll(7 + w/2, w). In 
the forward-jet problem, the energy scale is considered 
to be given symmetrically between the hard probes by 
log(W^VfcrQ) instead of \og{W^ / Q^) c::^ \og{l / x) as was 
the case for F2. In other words, we do not perform the 
shift 7^7-1-01/2 used for F2. 



Following Ref. , the renormalisation scale is fc^ ^ 
krQ- We have tested the sensitivity of our results when 
varying fc^ in the range [fc^-Q/A, A kxQ] with A = 2, with 
the substitution a{kTQ) ~^ a{XkTQ) + b a'^{kTQ) log{X). 
Note that, following formula the effective kernel is 
modified accordingly for each scheme, and we also modify 
the energy scale krQ^X kxQ- 

It is important to note that in formula we use the 
leading-order (Mclhn-transformed) impact factors given 
by 



-'^-2^^^2j^\ 27(1-7) 

q ^ 

r3(i + 7)r3(i-7) 

r(2-27)r(2 + 27)(3- 27) 



(5) 



as the full next-leading photon impact factors are not yet 
available (the jet impact factors are known at next-to- 
leading order [ij]). We point out that our phenomeno- 
logical approach can be adapted to full NLL accuracy, 
once the next-leading impact factors available. 

Expressing the integral in ([2]) using a saddle-point ap- 
proximation in 7, one finds for the theoretical forward-jet 
cross section 



a,(fc2)a,(Q2) 



feff{xj,k^)(l)l^^{-fc) 
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where x"// — ^X&f ! I 'H^ ■, and where the saddle point 
equation is Xgjj(7c,a) = 0. It is possible to extract the 
values of 7c, aXe//(7c,ci) and Q:Xe//(7c, after solving 
the implicit equation (|4]). They are given in Fig.l for the 
different schemes, as functions of a. 

The description of the forward-jet cross section is then 
almost parameter free; the value of a is imposed by the 
renormalisation group equations and only the overall nor- 
malisation is unknown (only the knowledge of the next- 
leading impact factors will provide a full prediction) . By 
comparison, the LL-BFKL formula is formally the same 
as (|ni), but with the substitutions 

Xeff Xll(7) = 2?/'(l) - V(l - 7) - "0(7) 
7c ^ 1/2 



a(fc2 



a. = const. 



with ■0(7) = '^logr(7)/c?7. In the LL-BFKL case, this is 
a two-parameter formula (normalisation and a) . The in- 
teresting property of our phenomenological approach is 
that formula ^ has formally the structure of the LL for- 
mula, but with only one free (normalisation) parameter 
and a NLL kernel. The delicate aspect of the problem 
comes from the scheme-dependent effective kernel Xeff- 
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FIG. 2: The forward-jet cross section da/dx measured by 
the HI coUaboration, compared with the three NLL-BFKL 
parametrizations S4, CCS and S3 using the krQ scale. For 
the S4 scheme, we also display the results of the 2kTQ and 
krQ/'i scales without changing the normalisation. 



FIG. 3: The forward-jet cross section da/dx measured by the 
HI collaboration, compared with LL and NLL (S4) BFKL 
fits and with NLO QCD predictions. The considered range of 
renormalisation scale for S4 does not alter the quality of the 
fit, after reajusting the normalisation. 



3. The NLL-BFKL formula for the fully differential 
forward-jet cross section is obtained from ([1]) and ([6]), 
with 7c, axe//(7c,<5) and ax"^y(7c,a) shown in Fig.l 
for the three different schemes. 

To fix the normalisation (the only free parameter) and 
check the quality of the data description using the BFKL 
formalism, we start by fitting the da/dx HI data 0. 
The choice of this data set corresponds to the kinemat- 
ical domain where the BFKL formalism is expected to 
hold (x <C 1 and Q^/k'^ ~ 1). The fitting procedure is 
the same as the one described in Ref. [1], Appendix A. 
The integrals over xj, and fcf. are performed numeri- 
cally taking into account the kinematic cuts given in Q. 
We performed fits on statistical errors only, systematics 
errors being strongly point-to-point correlated. In other 
words, we do as if the systematic errors were 100% cor- 
related (which is close to reality). 

The fit results to the da/dx HI data are shown in Fig. 2. 
The S4 fit can describe the data nicely (x^ = 5.4/5 d.o.f.) 
whereas the S3 and CCS schemes show higher values of 
^T^2 _ 45 5y'5 ^nd = 22.2/5 respectively). Note 
that, when fitting using total errors (shown on the fig- 
ures), all three schemes give values of below 1. In 
any case, the fitted normalisations keep reasonable val- 
ues, even with LL impact factors. We also display in 
Fig. 2 the results obtained when varying the renormalisa- 
tion scale in the range [kxQ/X^ A krQ] for the S4 scheme 
with A = 2. We notice that this change of scale essentially 
affects the overall normalisation and thus does not alter 
the quality of the fit, after reajusting the normalisation. 
For larger values of A, the quality of the fit deteriorates, 
due to the high values reached by a. 



A comparison of the S4 fit with the LL-BFKL results 
taken from Ref. @ is shown in Fig. 3. We notice the 
tiny difference between the LL and NLL results (the cor- 
responding curves are not distinguishable on the figure) . 
This confirms that the data are consistent with the BFKL 
enhancement towards small values of x. Contrary to the 
proton structure function F2, the forward-jet cross sec- 
tion da/dx does not show large NLL-BFKL corrections, 
once the overall normalisation fitted. This is due to the 
rather small value of the coupling a ~ 0.16 obtained in 
the LL-BFKL fit [Sj, corresponding to an unphysically 
large effective scale. 

We display in Fig. 3 the result obtained when the 7 
integration in ^ is computed exactly, and it is com- 
pared with the one obtained when using the saddle-point 
approximation ([5]). Considering the moderate values of 
y = 3 — 5 probed by the forward-jet measurement, one 
could question the validity of the saddle-point approx- 
imation. In the case of the S4 scheme, the comparison 
with the exact computation shows that, after all the kine- 
matic integrations have been performed, the difference 
can be absorbed in the overall normalization. This is also 
true for the S3 scheme. The case of the CCS scheme, for 
which there are difficulties to define the effective kernel 
in the whole complex plane, is left for future work. 

We also present in Fig. 3 the fixed order QCD calcula- 
tion based on the DGLAP evolution of parton densities. 
The next-to-leading order (NLO) prediction of forward- 
jet cross sections is obtained using the NLOJET-I--I- gen- 
erator fl5]. CTEQ6.1M [3| parton densities were used, 
and the renormalisation fij. and the factorization scale / 
were set equal to fj,^ = fij = Qk™""^, where fc™°^ corre- 
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FIG. 4: The forward-jet cross section da / dxdk^dQ^ measured 
by the HI collaboration, compared with LL and NLL (S4) 
BFKL fits and with NLO QCD predictions. The S4 resuR 
with exact integration of ([2} is also shown, and it generally 
agrees with the saddle-point calculation. 



sponds to the maximal transverse momentum of forward 
jets in the event. The NLO QCD predictions do not de- 
scribe the data at small values of x, as they are below by 
a factor of order 2. 

The fit parameters obtained with statistical error only 
are used in the following to make predictions for other 
observables. Namely, the relative normalisations between 
the different NLL-BFKL calculations (CCS, S3 and S4) 
are used to make predictions for the triple differential 
cross section da / dxdk^dQ'^ . This is an interesting observ- 
able as it has been measured with 3 different and 
cuts, yielding 9 different regions for the ratio r = k^/Q^. 
It was noticed in Q that the LL-BFKL formalism leads 
to a good description of the data when r is close to 1 and 
deviates from the data when r is further away from 1, as 
effects due to the ordering between Q and kr start to set 
in. NLO QCD predictions show the reverse trend. 

The HI data for da / dxdk^dQ^ are shown in Fig. 4 and 
compared with the S4 prediction, the LL-BFKL results 
(taken directly from [§]) and NLO QCD predictions. It is 
quite remarkable that the NLL-BFKL calculation, which 
includes some ordering between Q and kr, leads to a good 
description of the HI data on the full range. As it was the 
case for da/dx, the difference between the LL and NLL 
results is small when r ^ 1. By contrast when r differs 
from 1, the difference is significant, and the observable 
da / dxdk^dQ^ is quite sensitive to NLL-BFKL effects. 



When varying the renormalisation scale in the range 
IkrQ /2,2kTQ], the results shown in Fig. 4 differ by a 
small amount 17|, at most 10%. Also, we show in Fig. 4 
the curves obtained with the exact integration of and 
one can see that they are quite close to the saddle-point 
results, except for large values of r, where there is some 
deviation. Finally, we point out that the triple differen- 
tial cross section has the potential to resolve the scheme 
ambiguity, as the predictions of the other schemes CCS 
and S3 do not compare with the data as well as the pre- 
dictions of the S4 scheme 17 1. 

4. Let us summarize our main results. For the 
cross section da/dx, measured in the kinematical regime 
Q'^/k'^ ^ 1, the difference between the LL-BFKL and 
NLL-BFKL descriptions is very small (see Fig. 3), once 
the overall normalisation fitted. This confirms the va- 
lidity of the BFKL description of 0, Q previously ob- 
tained with the LL formula and a (rather small) effec- 
tive coupling. In the case of the triple differential cross 
section da / dxdk^dQ"^ , the same conclusion holds when 
r = k^/Q^ ^ 1. In addition when r differs from 1, the 
NLL-BFKL description is quite different from the LL- 
BFKL one, as it is closer to the NLO QCD calculation 
(see Fig.3). 

As a result, the best overall descrition of the data for 
da / dxdk^dQ^ is obtained with the NLL-BFKL formal- 
ism and the S4 scheme is favored. We need the com- 
plete knowledge of the next-leading impact factors be- 
fore drawing final conclusions, in particular with respect 
to the overall normalisations, but our analysis strongly 
suggests that the data show the BFKL enhancement at 
small values of x. This is of great interest in view of the 
LHC, where similar QCD dynamics will be tested with 
Mueller-Navelet jets [ll, [li|. 
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